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Section  1 


1.0  INTRODUCTION 

1 

In  a previous  publication,  we  have  shown  the  crucial  importance 

of  magnetic  fields  in  the  interpretation  of  continuum  x-ray  spectra 

16  2. 

obtained  when  high  Nd  laser  intensities  ( 10  W/cm  ) interact  with 
plane  solid  targets.  In  this  paper,  we  present  more  detailed  numerical 
results  from  a more  complete  model.  We  analyze  the  contributions  of 
the  various  thermal  magnetic  field  source  terms;  we  find  that  although 
the  magnetic  field  represents  only  a very  small  fraction  of  the  deposited 
laser  energy,  the  influence  of  the  B field  can  be  drastic  on  the  transport 
coefficients  and  can  in  turn  affect  the  overall  redistribution  of  energy 
in  the  plasma  produced  right  in  front  of  the  target.  Section  II 
describes  the  complete  B-field  equation  used  in  the  two-dimensional 
numerical  model.  Section  III  recalls  the  main  features  of  the  model 
and  the  initial  conditions;  Section  IV  consists  of  the  results  to  the 
B-field  equation  and  finally  Section  V summarizes  the  central  features 
of  the  solution  and  shows  its  influence  on  the  continuum  x-ray  spectrum. 


Note:  Manuscript  submitted  August  26,  1976. 
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Section  2 


2.0  MAGNETIC  FIELD  EQUATION 

The  B-field  is  derived  from  the  electric  field  which  is  itself 
derived  from  the  generalized  Ohm's  law.  All  terms  except  3 j/9 t are 
left  in  Ohm's  law  so  that  the  full  B-equation  reads: 

H * v*  [VxB  - fck'  <7xB>]-  % VTe 

' he  Vx  W~  (VxB)  xB  - | Vx  (i) 

e 

e 

where  V is  the  fluid  velocity,  ^ the  resistivity  tensor,  N0  the 

electron  density,  T the  electron  temperature  and  the  3 ' s are  coefficients 
e 2 

depending  on  u>ce  and  defined  in  Braginski  . The  only  assumptions 
made  in  deriving  Eq.  (1)  are  that  the  electron  pressure  is  isotropic  and 
that  the  ponderomotive  force  effects  have  not  been  included  in  the 
electric  field. 

The  last  term  on  the  right-hand  side  of  Eq.(l)  corresponds  to 
the  thermal  force  which  arises  due  to  the  unbalance  of  the  friction 
force  on  electrons  in  presence  of  a thermal  gradient.  This  term  enters 
the  expression  for  the  electric  field  as  - and  also  the  energy 

p /s  GING 

equation  as  - -j  where  = 3jVTe-B^  bxVTe. 

A short  analysis  of  the  terms  in  Eq.(l)  will  help  with  the 
interpretation  of  the  results  which  are  presented  further  down.  The 
third  term  on  the  right  hand  side  for  example  is  very  similar  in 
nature  to  the  second  one  if  we  neglect  the  term  due  to  curvature  of  the 
B-field.  It  can  then  be  written  as 

-in  w-'Jt) 


So,  only  when  B /8tt  becomes  large  compared  to  Pg,  this  term  becomes 
important. 


As  for  the  thermal  force  term,  because  the  B's  are  proportional 

to  the  B field,  it  can  be  seen  readily  that  the  spontaneous  magnetic 

field  will  grow  first  due  to  the  VN  x VT  term.  Only  when  (a)„,T  .) 

a e e ce  ei 

becomes  of  order  one,  the  thermal  force  terms  become  important  and  in 
fact  compete  with  the  electron  pressure  source  term.  It  is  expected 
then  that  the  thermal  force  term  alters  the  magnetic  field  significantly 
once  the  B field  has  reached  a magnitude  so  that  the  electron  gyrofrequency 
is  comparable  to  the  electron-ion  collision  frequency.  This  thermal 
force  term  being  a collisional  term,  the  condition  for  its  use  will 
be  that  the  electron-ion  collision  mean  free  path  must  be  smaller  than 
the  mesh  size.  If  we  take  for  0 


m • f • p 


1.44  10 


12  Te 


where  Tg  is  in  eV,  it  can  be  seen  that  \ will  be  less  than  12y  - the 


for  T <36  keV  at  the  critical 
e 


grid  size  used  in  the  calculation 
density  for  Z=6.  We  shall  see  that  this  condition  is  satisfied  there 
but  the  thermal  force  term  must  be  turned  off  in  the  underdense  region. 


Besides  the  magnetic  field  equation,  the  model  consists  of 
equations  for  the  conservation  of  mass,  momentum  and  energy  in 
cylindrical  coordinates  as  described  in  Reference  3.  The  target  is 
perpendicular  to  the  laser  axis  which  is  also  the  z axis.  Absorption  of 
laser  light  occurs  via  inverse  bremsstrahlung  for  electron  densities 
below  critical  and  at  the  critical  density  for  the  remainder  of  the 
laser  energy.  One  difference  with  Reference  3 is  that  the  present 
model  contains  a radiation  pressure  term  so  that  V*P  - e„  P„/c 
where  P{  is  the  absorbed  laser  power.  There  are  no  adjustable 
parameters  in  the  code. 

The  only  input  quantities  for  the  code  are  the  laser  pulse 
and  the  target  material  and  initial  configuration  at  t=0. 

The  laser  pulse  is  Gaussian  both  in  space  and  time  with  a FWHM 
of  21  psec,  an  energy  of  1J  and  a focal  spot  radius  of  12u.  The 
target  material  is  carbon  although  runs  with  polyethylene  were  made 
which  did  not  change  the  results  appreciably.  One  of  the  restrictions 
brouqht  by  the  use  of  an  eulerian  code  such  as  the  one  used  here  is  that 
the  calculations  cannot  start  with  an  infinitely  steep  gradient.  So  an 
initial  density  gradient  with  an  e-folding  length  of  27;i  was  set  up 
initially.  This  length  was  varied  from  20u  to  41p  without  bringing 
any  significant  changes  in  the  solution. 


Section  4 


4.0  RESULTS 

Profiles  for  the  magnetic  field  at  t=18,  36  and  88.5  psec 
are  shown  in  Figures  1 to  3.  The  thermal  force  is  not  included  in 
these  calculations  and  the  first  figure  corresponds  to  the  plasma 
variables  profiles  shown  in  Figure  3 of  Reference  1.  The  B field 
reaches  its  maximum  value  near  the  peak  of  the  pulse  - which  takes 
place  at  t=17  psec  - and  keeps  a value  close  to  its  maximum  until 
after  the  end  of  the  laser  pulse.  Changes  in  field  polarization 
occur  during  that  period  of  time  as  can  be  seen  from  the  figures  and 
are  due  to  intricate  structures  of  the  temperature  and  density  profiles 
in  two-dimensions.  It  is  also  seen  from  Figure  3 especially  that 
the  field  follows  the  fluid  motion  which  is  shown  typically  in  Figure  4 
and  5.  From  these  last  two  figures,  it  is  seen  that  the  plasma  expansion 
occurs  at  first  towards  the  laser  - which  is  located  at  the  bottom  of 
all  figures  - and  spreads  laterally  only  gradually  whereas  the  maximum 

O 

velocity  of  the  plasma  front  increases  to  speeds  of  10  cm/sec  (corresponding 
to  energies  of  62  keV  for  carbon  ions). 

In  Figures  6,  7,  and  8 the  magnetic  field  contours  corresponding 
to  the  same  case  are  shown  when  the  thermal  force  is  retained  in  Eq.  (1). 
From  the  comparison  between  this  series  of  figures  and  the  previous  ones, 
we  see  that 

a)  the  magnitude  of  the  maximum  field  is  reduced 
from  2.5  MG  to  about  1.3  MG. 

b)  the  maximum  B field  occurs  later  in  time  when  the 
thermal  force  is  included  and  at  a much  larger  radius 
(100u  vs.  12y ) 

c)  the  change  in  polarity  occurs  much  less  and  much 
later  for  the  thermal  force  case. 

Results  a)  and  c)  are  expected  from  our  analysis  of  the  effects 
of  the  thermal  force  terms,  namely  their  competition  with  the  electron 
pressure  source  terms  and  their  collisional  nature  leading  to  a more 
coherent  picture  overall. 


The  different  components  of  the  magnetic  field  source  terms 
can  be  written: 


3B 

fit 


c 3n  3T 
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and  are  shown  in  Figure  9.  They  correspond  to  the  case  of  Figure  6 
and  their  variations  are  shown  12p  off  axis, only  between 
z=320y.  Although  the  competition  between  them  is  strong. 


and  their  variations  are  shown  12y  off  axis, only  between  z=190  and 

3n  3T 

3R  3T  dz  3r 

the  dominant  term  followed  by  . 


is 


3r  3z 
3n  3T 
3r  -3z 

the  setting  of  the  artificial  density  gradient  in  front  of  the  target- 


Let  us  note  that  the 


term  which  comes  directly  from 


since  the  laser  steepens  the  density  profile  and  thus  produces  a 
positive  ~ - is  always  down  by  1 order  of  magnitude  from  the  dominant 
terms  and  in  that  respect,  the  initial  density  gradient  set-up  does  not 


invalidate  the  calculations.  Corresponding  to  these  different  various 

cases,  no  strong  maximum  temperature  off-axis  was  observed  as  has  been 
4 5 

reported  elsewhere  ’ . One  reason  for  this  result  is  that  in  our 

5 

own  previous  calculation  where  we  found  a maximum  off-axis  , the 
radial  grid  size  was  less  than  the  focal  spot  radius;  another  reason 
is  that  when  the  thermal  force  is  included,  the  maximum  B-field  is 
pushed  radially  outwards  and  heat  can  diffuse  more  easily  out  of  the 
focal  region. 


fi 


Section  5 


5.0 


CONCLUSION 


I ! 


In  both  cases,  with  or  without  thermal  force,  the  total 
magnetic  field  energy  is  several  orders  of  magnitude  (typically  six) 
below  the  absorbed  energy  and  it  is  interesting  to  note  that  although 
very  little  energy  is  converted  into  magnetic,  this  form  of  energy 
through  the  transport  coefficients,  control  the  repartition  of  absorbed 
energy  (thermal , radiated  and  even  kinetic  as  was  seen  in  Figures  4 
and  5). 

The  overall  effects  of  the  thermal  force  on  a macroscopically 
measurable  quantity  like  the  continuum  x-ray  spectrum  is  shown  in 
Figure  10  - along  with  the  corresponding  experimental  spectrum*.  The 
x-ray  spectrum  obtained  when  the  thermal  force  is  included  is  below 
that  one  obtained  without  the  thermal  force  and  then  agrees  better 
with  the  experiment.  This  improvement  can  be  explained  since  the 
reduction  in  heat  flow  is  less  when  the  thermal  force  is  included 
and  then  the  maximum  temperature  reached  in  this  run  was  19  keV  instead 
of  61  keV  for  the  corresponding  case  without  thermal  force.  The  dis- 
agreement above  100  keV  photon  energy  results  from  the  fact  that  the 
electron  should  be  treated  as  relativistic  when  it  reaches  that  energy 
level.  At  lower  photon  energy,  the  coarse  resolution  used  in  the 
calculations  did  not  bring  out  all  the  details  of  the  maximum  emission 
region  located  behind  the  critical  density  where  both  the  density  and 
temperature  gradients  are  very  steep  and  occur  over  one  cell  in  our 
calculations. 
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Fig.  3 — Still  later.  Note  complex  structure  and  multiple  field  reversals. 
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